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Abstract 



C^ ' The decay Kl — > n^^e'^e'^ occurs at a higher rate than the nonradiative process K^ -^ 

^ . vr'^e^e", and hence can be a background to CP violation studies using the latter reaction. 
It also has an interest in its own right in the context of chiral perturbation theory, through 

^vq . its relation to the decay Kl — > '/r'^77. The leading order chiral loop contribution to Kl -^ 

CN I 7r'^7e+e~, including the ((/e+ + qe-)'^/^^^ dependence, is completely calculable. We present 

^D ■ this result and also include the higher order modifications which are required in the analysis 

^ ; of Kl ^ vr°77. 
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I. Introduction 

There are three rare decay modes of the long hved kaon which have interrelated theoretical issues: 
Kl — > 7r°77, Kl — >• 7r°e+e~ and Kl — >■ 7r°7e"'"e~. The first two have been extensively studied; the 
latter has not been previously calculated. It is the purpose of this paper to provide a calculation 
of the latter process and describe how it is related to the phenomenology of the other two decays. 

There is a curious and important inverted hierarchy of these decay modes. The rate for 
the radiative decay Kl — »• 7r''7e"'"e~ is a power of a larger than the nonradiative transition 
Kl — > 7r°e~'"e~. This is because the Kl -^ 7r°e"'"e~ transition occurs only through a two-photon 
intermediate state, or alternatively through a one-photon exchange combined with CP violation 
(which numerically appears to be roughly of the same size as the two-photon contribution) |I|. 
The Kl — > n^e~^e~ rate is then of order a^. However, in Kl —>■ TT^'ye'^e~ we need only a 
one-photon exchange to the e~^e~ , leading to a rate of order a^. Our attention was first called 
to this inverted hierarchy by an observation that there are infrared divergences in a detailed 
study of the Kl — >■ 7r°e"'"e~ two-photon effect [jl| which need to be canceled by the one-loop 
corrections to the radiative mode Kl — > 7r°7e^e~ through the contributions of the soft radiative 
photons. This implies that the theoretical and experimental analyses of Kl — ^ ir^e'^e' and 
Kl — ^ 77r°e"'"e~are tied together. The soft and collinear photon regions of Kl -^ 'yir^e'^e' form 
potential backgrounds to the studies of CP violation in the Kl — > 7i^e~^e~ mode. 

The Kl — > 7r°7e"'"e~ mode also has an interest of its own. In recent years there have been 
important phenomenological studies of Kl —>■ 7r°77 in connection with chiral perturbation theory 
(ChPTh). This decay is calculable at one-loop {i.e., order E^*) ChPTh with no free parameters, 
yielding a very distinctive spectrum and a definite rate 0. Surprisingly, when the experiment 
was performed the spectrum was confirmed while the measured rate was more than a factor of 
2 larger than predicted. The way out of this problem appears to have been provided by Cohen, 
Ecker, and Pich (CEP) 0. By adding an adjustable new effect at order E^, as well as including 
known corrections to the Kl -^ vrvrvr vertex, they found that the predicted rate can be increased 
dramatically without modifying the shape of the spectrum much. This is also a surprising result, 
yet as far as we know it is the unique solution to the experimental puzzle. The ingredients of the 
mode studied in this paper, Kl — >■ 7r°7e"'"e~, are the same as for Kl — >■ 71^77, except that one 
of the photons is off shell. Within the framework of the CEP calculation, the ingredients enter 
with different relative weights for off-shell photons. This will allow us to test the consistency of 
the theoretical resolution proposed for Kl — »■ 7r°77. 

We outline the computation for the (9(E'^) contribution to the process in Sec. II, and then 
we extend it to C(E^) in Sec. III. Finally, we recapitulate our conclusions in Sec. IV. 

II. The OCE^) calculation 

First let us provide the straightforward 0{E^) calculation within ChPTh. This is the general- 
ization to kf y^ of the original chiral calculation of EPR [Q . Here ki is the momentum of the 
off-shell photon. This captures all the k\/rn^ and k\/m\ variations of the amplitudes at this 
order in the energy expansion. There can be further k'l/{l GeV)^ corrections which correspond 
to (9(E^) and higher. The easiest technique for this calculation uses the basis where the kaon 
and pion fields are transformed so that the propagators have no off-diagonal terms, as described 
in ref. [^. The relevant diagrams are then shown in Fig. 1. Defining ^ as 

g = Gg/S, Gg = GplV^dVuslgs, \9s\ ~ 5.1, (1) 

the diagrams give the following integrals, respectively: 
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Figure 1: Diagrams relevant to the process Kl -^ 7r°7e^e^ at 0{E'^) and 0{E^) 
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Interestingly when we add these together the i^ ^ Svr amplitude factors out from the remaining 
loop integral resulting in 
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It is not hard to verify that this result satisfies the constraints of gauge invariance k^M.^^ = 
k2-Mfj_i, = 0. At this stage, the integral may be parametrized and integrated using standard 
Feynman-diagram techniques. Let us keep photon number one as the off-shell photon and set 
/c| = 0. In this case the amplitude with one photon off-shell is described by 
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The notation is defined by 



s= {PK-Pof = {ki + k2f 



and 



^ , , f^ dzi ^ 
F(a)= / ^log 
Jo Zi 



G{a) = / dzi log 
Jo 



9 

mi 


-a{l-zi)zi 


— ie 




ml 




ml 


-a{l - zi)zi 


— ie 



mt 



The above functions are related to those presented by CEP 0]: 
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This agrees with the EPR result in the kf ^ limit. 
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Figure 2: Differential branching ratio dr{KL -^ 7r°7e"'"e )/dz to order (9(E^). 

At this order we have also calculated the additional contribution resulting from the kaons 
circulating in the loops of Fig. 1. They give rise to 
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The resulting integral is similar to that of Eq. (8), substituting the mass of the pion with that 
of the kaon. Attaching an e^e~ couple to either photon and adding all the above contributions 
together, the result we obtain for the branching ratio is 
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the decay distributions in z and y provide more detailed information. We present them in Figs. 
2 and 3. 



III. The 0{E^) calculation 



We also wish to extend this calculation along the hues proposed by CEP 0, who provide a 
plausible solution to the problem raised by the experimental rate not agreeing with the (9(E^) 



dBR(KL^7T V e")/dy 



3.010 




Figure 3: Differential brancliing ratio dr{KL —>■ 7r°7e"'"e )/dy to order (9(E^). 

calculation when both photons are on-shell. The two primary new ingredients involve known 
physics which surfaces at the next order in the energy expansion. The first involves the known 
quadratic energy variation of the K ^ Sn amplitude, which occurs from higher order terms 
in the weak nonleptonic Lagrangian [^, §. While the full one-loop structure of this is known 
0, it involves complicated nonanalytic functions and we approximate the result at (9(E^) by 
an analytic polynomial which provides a good description of the data throughout the physical 
region: 
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ai and 02 are obtained from a fit to the amplitude for Kl -^ vr tt+tt^ [Q and to the amplitude 
and spectrum for Kl -^ vr^e^e" 0], so that their values are constrained within their theoretical 
uncertainty of 10 - 20%. We have numerically verified that such a variation of said parameters 
involves a very modest change in the shape of the spectrum for K^ —>■ 77r°e"'"e~ and a change in 
its final branching ratio somewhat smaller than the uncertainty on the parameters. 
The other ingredient involves vector meson exchange such as in Fig. 4. Some of such contributions 
are known, but there are others such as those depicted in Fig. 5 which have the same structure 
but an unknown strength, leaving the total result unknown. In Ref. the result is parametrized 
by a "subtraction constant" which must be fit to the data. 

In principle one can add the ingredients to the amplitudes and perform a dispersive calculation 
of the total transition matrix element. In practice it is simple to convert the problem to an 
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Figure 4: Vector meson exchange diagrams contributing to Kl — * 7r°7e"'"e~. 

effective field theory and and do a Feynman diagram calculation which will yield the same 
result. We follow this latter course. 

The Feynman diagrams are the same as shown in Fig. 1, although the vertices are modified 
by the presence of (9(E^) terms in the energy expansion. Not only does the direct K ^ 3n vertex 
change to the form given in Eq. (18), but also the weak vertices with one and two photons have 
a related change. The easiest way to determine these is to write a gauge invariant effective 
Lagrangian with coefficients adjusted to reproduce Eq. (18). We find 
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The resulting calculation follows the same steps as described above, but is more involved and is 
not easy to present in a simple form. We have checked that our result is gauge invariant and 
reduces to that of CEP in the limit of on-shell photons. 

The contribution proportional to ai can be computed analogously to those already calculated 
for the O(E^) case: 
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Figure 5: Vector meson exchange diagrams contributing to Kl 
strength. 
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From the above formulas we obtain 



M 



liv 



(47r)2 

+ B{xi,X2 



A{xi,X2){k2^kiy - ki ■ k2g^,v) 
Pk ■ hPK ■ ^2 



ki ■ k2 



+ D{x„X2){kf-f^ 
\ ki-k2 



9nv 



-kifj,kiu + kifj_pj^j^ — - ^k2^Px,/ 



ki ■ k2 



ki ■ k2 



where 



(27) 



(2^ 



A = 16a2e^{2[l - 2(a;i + X2)]Ii{ziZ2) + Xili{z2) + X2[2Ii{zl 

- 32a2e'{[2xl - xi{z + q)][-l2izfz2) + h{zlz2)] 

+ [2a;ia;2 - Xi{z - q)/2 - X2{z + q)/2][2l2{zlz^) + l2{ziZ2) - 12(2? ^2 

- hizizl)] + [2x2 _ ^^(^ _ g)][/2(^,4) _ l2{z,zl)]} 



hiz2)+lM)]} 



A 777 

+ -a2e2log^ + (47r)2VMD^, 
3 mj 



(29) 



4 / 771^ \ 

B = -32a2e'/3 + 1602/4 + -a2e\z - g) -1 + log ^ + {AnfVMDB, 

6 \ mp) 



(30) 



with 



D 



- I + 16a2e2[2a;2 - {z - q)/2][2h{ziZ2) - I^)] 

+ lQa2e\2y-q)[h{z^)-h{l)/2]+Aa2e^[2x^ - {z + q)/2]h 
+ (47r)2VMD,5, 



h{z-z^) 



1 rl-zi £) 

dZi / C?Z2-2"-2^l0g^, 



JO 



mt 



(31) 



(32) 



H^lZT) = dzJ dZ2^ 

JO Jo Ui 



h = dzi dz2Di log -^ 

Jo Jo mi 



h= dziDilog—^, 
Jo m^ 

h= [ dzi{Azl ~ Azi + l)\og 
Jo 



mi 



and 



D2 

Xi 



ml - 2ki ■ k2ZiZ2 - klziil - Zi), 

ml -klzi{l- zi), 

Pk -ki Pk- h 



m 



X2 



K 



m 



K 



VMDa(Xi,X2) 



V=iij,p 



Pk ■ (Pk - h) ^ Pk " (Pk - K 



{Pk-^^Y -"^v {Pk-^iY -'^v 



V=LO,p 



+ 



{Pk - ^2)^ - mlr {pk - hy - ml 



VMDd{xi,X2)= J2 Gv 
assuming the numerical values 



ki ■ ^2 



v=u,p {PK-kiY-my 



Gpm\ = 0.68 X 10 



-8 



G^m\ 



-0.28 X 10 



-7 



(33) 
(34) 
(35) 
(36) 



(37) 

(38) 
(39) 
(40) 

(41) 



The loop calculation that we have just described provides all of the off-shell dependence 
scaled by the pion mass, and is of the form kf/ml. There can be an additional dependence 
of the form k\/ A?, where A ?^ 1 GeV. We cannot provide a model independent analysis of the 
latter. However, experience has shown that most of the higher order momentum dependence is 
well accounted for by vector meson exchange. Therefore we include the k\/ A? dependence which 
is predicted by the diagrams of Fig. 4. One can recover the parametrization in ay neglecting the 
dependence on {pK~kiY and (pj^ — /c2)^ in formulas (38) - (40), and performing the replacement 
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where G^s ~ Gp+G^. This completes our treatment of the Kl -^ 7r°7e"'"e amplitude. 
The calculation we have presented in this section leads to the total branching ratio of 

BR{Kl ^ 7r°7e+e-) = 2.3 x 10"^ 
The decay distributions are presented in Figs. 6 and 7. 
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Figure 6: Differential branching ratio dT{KL -^ 7r°7e"'"e )/dz to order (9(E^). 

IV. Conclusions 



7r*^77. 



The behavior of the Kl -^ tt 7e+e~ amphtude mirrors closely that of the process Kl 
The more complete calculation at order E^ gives a rate which is more than twice as large as the 
one obtained at order E*^, despite the fact that the new parameter introduced at order E^ is quite 
reasonable in magnitude. This large change occurs partially because the order E^ calculation is 
purely a loop effect, while at order E^ we have tree level contributions, and loop contributions 
are generally smaller than tree effects at a given order. It was more surprising that the spectrum 
in Kl -^ 7r°77 was not significantly modified by the order E^ contributions. These new effects 
are more visible in the low z region of the process we have calculated, K^ — * n^'-fe~^e~ . 

This reaction should be reasonably amenable to experimental investigation in the future. It 
is 3-4 orders of magnitude larger than the reaction Kl -^ Ti^e^e~ which is one of the targets of 
experimental kaon decay programs, due to the connections of the latter reaction to CP studies. In 
fact, the radiative process of this paper will need to be studied carefully before the nonradiative 
reaction can be isolated. The regions of the distributions where the experiment misses the photon 
of the radiative process can potentially be confused with Kl — > 7r°e"'"e~ if the resolution is not 
sufficiently precise. In addition, since the 7r° is detected through its decay to two photons, there is 
potential confusion related to misidentifying photons. The study of the reaction Kl -^ 7r°7e"'"e~ 
will be a valuable preliminary to the ultimate CP tests. 
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Figure 7: Differential brancliing ratio dr{KL —>■ 7r°7e"'"e )/(iy to order 0{E^). 
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